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Images of the Hubble Ultra Deep Field are analyzed to obtain a catalog of galaxies 
for which the angular sizes, surface brightness, photometric redshifts, and absolute 
magnitudes are found. The catalog contains a total of about 4000 galaxies identi- 
fied at a high signal-to-noise ratio, which allows the cosmological relations angular 
size-redshift and surface brightness-redshift to be analyzed. The parameters of the 
evolution of linear sizes and surface brightness of distant galaxies in the redshift in- 
terval 0.5-6.5 are estimated in terms of a grid of cosmological models with different 
density parameters (r^y; ^Im)- The distribution of photometric redshifts of galax- 
ies is analyzed and possible superlarge inhomogeneities in the radial distribution of 
galaxies are found with scale lengths as large as 2000 Mpc. 



I. INTRODUCTION 

The program of observational cosmology, which was first formulated by Hubble and 
Tolman [ij and then further developed by Sandage [21, 0] suggested a number of cosmologi- 
cal tests — N(m), N(z), m{z), J{z), Q{z), and t{z) — using numbers of objects, magnitudes, 
surface brightness, sizes, and their ages. These tests, which are also called classical tests, 
are based on the comparison of empirical relations between directly observable quanti- 
ties with the theoretical relations between the same quantities as predicted by different 
cosmological models. 

Modern approach toward the analysis of classical cosmological tests consists in the 
simultaneous taking into account both the parameters of the cosmological model and the 
evolution of galaxies. However, so far, no bona fide model of the evolution of galaxies 
is available and the development of such a model remains the main unsolved problem of 
modern cosmology. 

The evolution of galaxy parameters is usually estimated for the "standard" values of 
cosmological parameters flm = 0.3 and fly = 0.7 with reference to WMAP data. However, 
as Spergel et al. [3] pointed out, interpretation of observations of the fluctuations of 
microwave background includes 15 parameters of the standard model and only six of 
them can be estimated independently. Note that the parameter Qy of vacuum density 
is not among those six parameters. Parameter Qy is estimated using a combination of 
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other observational data, such as the Hubble diagram for type la SNe and the correlation 
properties of the large-scale distribution of galaxies. 

Note that the values of cosmological parameters estimated using different methods 
may differ significantly from the corresponding "standard" values. Indeed Clochiatti 
et al. used SNIa observations in the "Hight-Z Supernova Search" program to es- 
timate the following parameters: = 0.79 ± 0.15 and Qy = 1.57 ± 0.25. The re- 
sults obtained in "The Supernova Cosmology Project" [sl using the light curves of 
type la supernovas obtained simultaneously in several filters yielded the estimates 
Qm = 1-26 ± 0.4, Qy = 2.20 ± 0.5. Thus the estimates of cosmological parameters Qm 
and Qy may vary over a wide range. 

In this paper we perform a quantitative analysis of the dependence of the parameters of 
galaxy evolution on the cosmological model parameters Qm and Qy We use results of ob- 
servations of the Hubble Ultra Deep Field to obtain the diagrams "angular size-redshift" 
{Q{z)) and "surface brightness-redshift" {J{z)) for galaxies in the redshift interval 0.5- 
6.5 and analyze the effect of a change of cosmological parameters on the estimates of the 
galaxy evolution parameters. 

In addition, we also analyze the distribution of photometric redshifts of galaxies — the 
dN{z)/dz test — for a magnitude-limited sample. Our method is capable of identifying 
large-scale fluctuations of the number of galaxies, which exceed the Poisson noise level. 



II. OBSERVATIONAL DATA 

Problems involving cosmological tests require reliably determinable galaxy parame- 
ters, i.e., to address such problems, one must use values of the numerous parameters of 
photometric reduction and identification of distant galaxies and ensure sufficiently high 
signal-to-noise ratio. To do this, we compiled a catalog of galaxies of the Hubble Ultra 
Deep Field, including objects with the signal-to-noise ratio higher than 5. 



A. Identification of Galaxies 

The data on the Hubble Ultra Deep Field (HUDF) can be found at 



http://www.stsci.edu/hst/udf. We used the images taken in four filters adopted from 



http://archive.stsci.edu/prepds/udf/udfMsp.html, where the initial reductions have already 



been performed, allowing the identification of objects on frames to be immediately ad- 
dressed. We used sExtractor [3] code to identify objects. Our input data consisted of 
a configuration file and files with the frames taken in four filters (B, V, i, and z). In 
the configuration file we set the parameters that we used to identify objects in the field 
studied. 

We set the PIXEL_SCALE parameter equal to 0.03" for the field considered. The pro- 
cess of identification of objects depends significantly on parameter DETECT_THRESH. 
sExtractor interprets the signal as a part of the galaxy if the flux level exceed this param- 
eter. 



Our criteria of the identification of objects are based on the following assumptions: 



1) all pixels record signals that exceed the given DETECT_THRESH; 

2) pixels form a group (i.e., they are "crowded"); 

3) the number of pixels in a given group is greater than the given natural number. 

We considered the pixels with flux deviated by more than 3a from the mean flux to be 

parts of the object. 

The "crowdnectess" of pixels was set by parameter DETECT_M1NAREA. If a group 
of pixels has a count above 3a and the number of pixels exceeds DETECT_MINAREA, 
the programm concludes that a galaxy has been detected. 

The mean value is estimated by averaging a portion of the fleld rather than the entire 
field. The given averaging region had the size of 100 x 100 pixels. The finding of the 
background averaging region is an issue of great importance, because use of too small 
regions effects an estimate being obtained, and use of large regions results in a strong 
effect on the detected objects. 

We also additionally smoothed the image. The smoothing procedure is performed 
before the detection of galaxies of the field. In this paper we use Gaussian as the smoothing 
function. The size of the smoothing region is 3 x 3 pixels. In particular, the 3x3 
region has the form of a square matrix (the axes X and Y), where each cell stores the 
corresponding normalized flux (the Z coordinate). In this way the three-dimensional 
model of the smoothing function is modeled, where X and Y correspond to the coordinates 
in pixels and Z, to the flux. 

In cases of the search and separation of objects we used an additional parameter, which 
allows to take into account the input from the bright objects. All detected objects are 
tested for the closeness to bright objects. In cases of "contamination" by a bright object 
we use the following formula of the Moffat proflle: 

m = I m 

J(0) {l + kxr^y ^ ' 

where /3 is set in the sExtractor conflguration flle and J is the surface brightness. A 
change in parameter /3 affects signiflcantly the photometry and detection of fleld objects. 

B. Finding the Parameters of Identified Objects 

We found the main parameters of the objects identifled. 

1. Photometric •parameters 

The photometric parameters are characterized by the flux from the given area of the 
object. There are several possible methods of identifying the area the flux is to be com- 
puted from. In this paper we use the "isophotal" approximation set in the sExtractor 



configuration file. In this method the contours of the area are found from the count level 
that depends on the average flux in the given area. Figure [T] shows the isophote of one 
of the galaxies in the field studied. The background is subtracted from the flux gathered 
from the given area and the magnitude of the detected object is computed. 
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FIG. 1: Isophote of one of HUDF galaxies. 

While compiling the catalog of objects we tried to use different methods of flux com- 
putation and found no significant differences between the results of the photometry of 
galaxies (depending on the method). This may be due to the fact that most of the 
galaxies in the field studied have small sizes and irregular structure. 

The instrumental flux is related to magnitudes via the following formula: 

m = 2.5 X lg(F) +m^p, (2) 

where F is the flux in instrumental units and rrizp, the average magnitude of background 
for each filter. We also find the maximum surface brightness of the object, which is 
available via parameter MU_MAX in mag" 

The input catalog of objects contains the following photometric information: 

• The instrumental flux with its error; 

• The apparent magnitude of the object with its error; 

• The maximum surface brightness of the galaxy; 

• The effective radii corresponding to 25%, 50%, and 75% of the flux of the entire 
galaxy. 

2. Astronomical parameters of galaxies 

We use the barycentric coordinates of the objects computed by sExtractor in accor- 
dance with the following formula: 



X = .= ^^, (3) 

where /j are the moments equal to the galaxy flux in each pixel. 

The coordinates of the objects are computed both in the equatorial coordinate system 
for the epoch of 2000.0 and in the relative (Cartesian) coordinate system. The relative 
coordinate system is determined by the coordinates of galaxies in terms of image pixels. 



3. Geometric parameters 

The geometric parameters describe the sizes and appearance of the objects. The el- 
lipticity of galaxies is characterized by the semiminor and semimajor axes (o and b, re- 
spectively) and by position angle Qse- The semiminor and semimajor axes are computed 
using second-order moments. 

The formulas for the second-order moments have the following form: 
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The semiaxes can then be computed by the following formulas: 
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The position angle Qge is counted from the North direction for the epoch of 2000.0. 
The oblongness of the object is characterized by its elhpticity (1 — |) or elongation (|). 
We also compute the area of the object at the level of the detection threshold DE- 
TECT_THRESH. 



4. Sizes of galaxies 



Cosmological tests usually employ the effective sizes of galaxies at half the level of 
flux proflle — the so-called FWHM (Full Width at Half- Maximum). We described the flux 
proflle by a Gaussian. 



C. Identification of Objects 



After setting the parameters in the sExtractor configuration file we created a file used 
to identify objects in the entire field. We obtained it by coadding the FITS files (frames 
of the field studied taken in four filters) using MIDAS package. We did it because objects 
appear greater in the i and z filters compared to their sizes in the B and V filters. To 
compile the preliminary catalog in each filter, we use the file of object detection and the 
file with the image. As a result, we found more than 4300 objects. 

We then produced a sample of objects by imposing the following constraints. First, 
the signal-to-noise ratio (S/N) must be no lower than 5. We excluded from the sample 
all objects with S/N smaller than 5. We further excluded all objects with no measured 
fiux in at least one of the four filters. As a result, we obtained a catalog containing 4125 
galaxies. 

D. Finding the Redshifts 

After identifying the objects and compiling the preliminary catalog in four filters we 
found the photometric redshifts z and absolute magnitudes of galaxies. 

The photometric redshifts of galaxies are inferred from their magnitudes in different 
filters. The efficiency of this method is based on identifying the photometric data points to 
portions of the continuum spectrum of the galaxy. The accuracy of photometric estimates 
is lower than that of spectroscopic estimates and depends on the set of filters employed 
and the accuracy of photometric data. However, photometric redshifts are quite suitable 
for many cosmological and extragalactic problems. This method of redshift finding has 
been playing important part in observational cosmology. 

Although finding of photometric redshifts requires no spectroscopic data, we must, 
however, trust the template galaxy spectra used for comparison. 

We use HyperZ code jsf to compute the photometric redshifts. 

The input parameters of HyperZ code include: 

• Apparent magnitudes of objects in the filters B, V, i, and z and their errors; 

• Template spectral energy distributions (SED); 

• The reddening law for the objects; 

• Cosmological parameters. 

1. The method of redshift finding 

The procedure of the finding of photometric redshifts is based on the comparison 
of the observed energy distribution in the spectrum of the galaxy with the template 
spectral energy distribution. The observed distribution (inferred from photometric data) 
is compared to different template galaxy spectra using the same photometric system. 



We find the photometric redshift z for the given object from the template spectrum 
that fits best the observed spectrum. The fitting procedure is based on minimization of 
X^. We compare the observed and template distributions by the following formula: 

2 1 \ [FobsA b X Ftemp,ii,Z^ l'^ /„\ 

X {z) = 2^ , (9) 

i=l 

where Fobs,i is the observed flux; Ftemp,, the template flux; (Xj, the variation of the flux in 
the given filter, and b, the normalization constant. We thus find the redshift by minimizing 

HyperZ code includes both observed |9[ and synthetic template spectra. Template 
spectra can be replaced and are stored in the configuration file. 

As the synthetic models of galaxy spectra we use GISSEL98 datasets (Galaxy Isochrone 
Synthesis Evolution Library), which feature spectra from a wide interval of energies: from 
UV to far IR with the possibility of evolution to large z. 

Note that we computed the photometric redshifts using template redshifts for galaxies 
of the E/SO, Sa, Sb, Sc, Sd, Im, and Burst types. 



2. Various corrections and reductions 

Metallicity. HyperZ code can take into account the evolution of galaxy metallicity. 
To do this, special template galaxy spectra must be indicated that take this factor into 
account. We found the photometric z using template spectra with and without the ac- 
counting for metallicity evolution. The results obtained led us to conclude that metallicity 
evolution has no effect on the inferred photometric redshifts. That is why we did not take 
metallicity evolution into account while finding z. 

Reddening law. Recent studies of galaxies at large z demonstrated the importance of 
the accounting for the "reddening" when estimating the redshifts because of the effect 



of Galactic dust in the object studied. The Calzetti law [10|] describes fairly well the 



extinction for objects at large redshifts and we adopted it in this paper. Calzetti et 



al. 10|] describe the empirical results for the group of young galaxies, which we used for 



finding z: 



k{X) = 2.659 f-2.156 + - + + Ry, 0.12/im < A < 0.63/xm 

k{X) = 2.659 (^-1.857 + + Rv, 0.63fim < A < 2.20fim , (10) 

where Ry = 4.05 — is the total extinction in the V-filter. 

This law corresponds to central star-forming regions and it can therefore be applied 
to galaxies at large redshifts. Note that extinction is important in the UV part of the 
spectrum. As a result, extinction must show up appreciably for galaxies with z > 3 in 
the optical part of the spectrum where the UV portion has been shifted. 



Further corrections. The radiation of distant galaxies is also "distorted" by extinction 
in HI regions located along the line of sight, which shows up in the part of the spectrum at 
wavelengths shorter than Lya(A = 1216A). HyperZ code allows for this effect and applies 
corrections in accordance with the law suggested by Oke and Korycansy 11 



3. Parameters used to find the redshifts 



To find the photometric redshifts, HyperZ code needs a configuration file with the 
following data: 

• Different cosmological models with Hq = 72 km/s x Mpc~^ for computing absolute 
magnitudes M in the B filter; 

• The magnitudes from the input catalog of galaxies in various filters; 

• The transmission parameters of the B, V, i, and z-band filters in accordance with 
HST data; 

• The redshift step, Az = 0.1. 



E. Finding the Surface Brightness Profiles of the Galaxies 

The surface brightness profiles of galaxies are computed via iterations by the following 
formula: 



X 



2 ^{fl^-^X,y moddx^y) ' (11) 



x=l y=l '^x,y 



where fluxx,y and modelx,y correspond to the observed and model flux, respectively and 
cr to the background weight. We minimize \f over the entire image (the coordinates x 



and y) and use GALFIT for computations 12 



We fit the surface-brightness distribution by a Gaussian. It would make no sense to use 
a more complex model for the surface-brightness distribution, because most of the galaxies 
in the field studied do not fit standard classification and have irregular appearance. 

We adopt the model surface-brightness profile in the following form: 



.2 



= Jo X exp (^) , (12) 

where Jq is the central surface brightness; r, the galactocentric distance, and a = 
FWHMsb/2M5. (here FWHMsb is the full width of the Gaussian at half maximum). 
We compute Jq by the following formula: 

Ftot X R{c) 

- 2na^q ' ^^^^ 



where Ftot is the total flux from the object studied and q, one of the input approximation 
parameters of GALFIT code, 

X 7r(c + 2) 

''^'^-4/?(l/(c + 2),l + l/(c + 2))' 

(3 is the Beta function, and c is the galaxy ellipticity parameter. 

We set the following input parameters for finding of the surface-brightness profiles: 

• The size of the field containing the object (a FITS file); 

• The size of the field to be fitted; 

• The background level in magnitudes; 

• The a value of the background for each pixel in the given region of the field (the 
weight image, a FITS file); 

• The pixel scale in arcseconds. 



The initial data (the zero approximation of the theoretical model of surface-brightness 
distribution) include: 

• The coordinates of the center of the object; 

• The total fiux from the entire object in magnitudes; 

• The FWHM of the surface-brightness profile of the object; 

• The axial ratio b/a of the galaxy; 

• The position angle of the galaxy. 

We computed the surface-brightness profile only for galaxies with absolute magnitudes 
M ranging from -20 to -18 and with the 99-100% probability of photometric rcdshifts 
(zphot) (it is one of the parameters in the general catalog of objects). After selecting 
objects from the catalog in accordance with the above criterion we studied the objects 
with FWHM flux > 10 pixels. We imposed this additional constraint on the sample in 
order to improve the convergence of the method that we use to find the surface-brightness 
profile. 

As a result, wc obtained a catalog of parameters for finding the surface-brightness 
profiles in the Gaussian approximation for each object of the sample considered. 
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FIG. 2: The distribution of the observed quantities of the galaxies of the catalog for the 
CiJbviz-modS model, where (a) is the distribution of photometric redshifts; (b) — the distri- 
bution of B-band absolute magnitudes; (c) — the distribution of i-band apparent magnitudes, 
and (d) — the distribution of angular sizes. 

III. ANALYSIS OF THE GALAXY SAMPLE 
A. Constructing the Main Catalog 

We already mentioned above that after reducing the HUDF we constructed the general 
catalog of the objects. However, we then modified this catalog and drew various samples 
from it. To systematize the data, here wc review the samples and catalogs. 

After the objects were identified by sExtractor, we constructed the preliminary catalog 
of data, CiJbviz. It contains the principal parameters (the coordinates, sizes, ellipticity, 
fluxes, etc.) for each object in the four fllters B, V, z, i. 

To flnd the photometric redshifts, we compiled the input catalog of data containing 
the magnitudes of the galaxies in the four fllters and their color indices. We added the 




FIG. 3: The "absolute magnitude-linear size" relation for all galaxies of the sample CiJbviz-modS 
with the parameters of standard cosmological model. 

input data including Zph (photometric redshifts), M (absolute magnitudes), and Galaxy 
Types into the the general data catalogs CiJwizjmodl, CiJt)vizjmod2, CiJbvizjmodS, and 
CiJbvizjmodA corresponding to the cosmological models, — 0.0 — 0.0, = 0.0 
Qjn — 1-0, — 0.7 — 0.3, and Qy — 1-0 — 0.0. We then subscripted the name 
of the catalog depending on the absolute magnitude of objects selected from the main 
catalog: M18_20 (galaxies with absolute magnitudes M in the interval from -20 to -18) 
or M20_22. 

We the excluded from the resulting catalog of objects the galaxies with incorrectly 
computed flux in at least one of the fllters (in these cases sExtractor outputs the the 
flux value equal to 99 magnitudes) and the probability of flnding the photometric redshift 
Zph < 90%. 



B. Distribution of the Observed Parameters of the HUDF Galaxies 



In Fig. [2] we represent the distribution of the observed parameters of the objects of the 
i-band catalog. The histograms of 0, m, M, and z are computed for the galaxy sample 
of the CiJ)viz_mod3 catalog, which contains a total of 4125 objects. 

It is evident from Fig. [2] that the distribution of photometric redshifts is significantly 
inhomogeneous (a); the distribution of absolute magnitudes M lies in the interval from 
-24 to -12 (b); the number of galaxies increases with decreasing flux down to 29"^, which 
corresponds to the completeness limit of the catalog (c); the characteristic angular size B 
is equal to 0.3" (d). 

C. Constraints on Absolute Magnitudes 

Figure [3] shows the "absolute magnitude — linear size" diagram for galaxies. We use 
the FWHMfiux quantities, which correspond to the halfwidth of the flux profile, as the 
sizes of the corresponding galaxies. 

The increase of linear size becomes slower for bright galaxies and therefore, to reduce 
the systematic displacement of the data points on the Q{z) diagram, constraints are to 
be imposed on the absolute magnitudes of galaxies. The distribution of linear sizes of 
faint galaxies is uniform and therefore we exclude them in order to reduce the system- 
atic errors. The vertical line in FigjH] corresponds to the limiting absolute magnitude 
of galaxies included into the sample, which contains galaxies with absolute magnitudes 
brighter than -18. We compute the absolute magnitudes using formula (13T!) . which al- 
ready incorporates a fixed cosmological model and spectral energy distributions of galax- 
ies. Therefore to estimate the evolution parameter k, we use the subsamples of galaxies 
limited by absolute magnitude in accordance with each particular cosmological model. 

D. Division into Subsamples 

We found the parameter of the evolution of linear sizes of galaxies as a function of 
their type and subdivided the sample into the subsamples of spiral and elliptical galaxies. 
HyperZ code determine the galaxy type by the form of its spectrum while computing 
the photometric redshifts. We selected the brightest galaxies of the catalog for various 
intervals of absolute magnitudes. We computed the absolute magnitude M separately for 
each cosmological model. Then we subdivided the CiJwiz catalog into four subsamples 
(in accordance to each cosmological model), which we then subdivided into subcatalogs: 
-20 < M < -18; -22 < M < -20. 

IV. THEORETICAL e{z) AND J{z) RELATIONS 

Many papers have been dedicated to the discussion of observational tests of cosmo- 
logical models (see, e.g., the reviews by Sandage [I, lil). Modern data are indicative of 



the need for the use of cosmological models including both dark matter and dark en- 
ergy. In this section we give the general theoretical relations used to analyze the classical 
cosmological tests. 



A. The "Metric Distance Redshift" Relation 



The density parameter in the Friedmann cosmological models including cold dark mat- 
ter (CDM) and dark energy is equal to the sum 

n = nm + ^v, (15) 

where flm = Pm/ Pcru is the parameter of the density of cold matter and fly = Pv/Pcrit is 
the parameter of the density of dark energy, for which pv = wpvc^, and w < 0. In the 
particular case w = —1 dark energy corresponds to cosmological vacuum and Einstein's 
cosmological constant. 

In the general case the Friedmann equation has the following form: 

Q = l-Qk, (16) 



or 

2 SttG kc'^ 
H -^P = -^. (17) 

where Q = g/ Qcrit is determined by the total density Q = Qm + Qv] the critical density is 
Qcrit = 3if^/87rG, and the curvature density parameter is flk = kc^/S'^H'^. The Hubble 
parameter H = S/S and the scale factor S{t) are set by the Robertson- Walker four- 
interval, which has the following form: 

ds^ = c^de - S\t)dx^ - S\t)ll{x){de^ + sin^ ed(P'') , (18) 

where Ikix) = sin(x), x, sinh(x) with k = +1, 0, and -1, respectively. 

The proper metric distance r from the observer to the galaxy with dimensionless co- 
moving coordinate x metrics (ITSi) is given by the following formula: 

r{t,x) = Sit)x. (19) 

Note that the proper metric distance r (measured inside the three- 
dimensional hypersphere) and scale factor S have the dimensions of length: 

[r] = [S] = [cm] . 

To describe the "angular size-redshift" relation, observational cosmology uses the 
"external" metric distance / 

l{t,i^) = Sit)fx, (20) 
where dimensionless comoving distance /i appears explicitly in the four-interval 

ds' = c'dt' - S\t)-^^ - S{tfp\de^ + sin^ edcp"). (21) 

1 /b^ 



According to formula (I2T1) . distance / = Sfi is measured in the ambient four- dimensional 
space. The relation between x and fi {x = ^k^ilAi = -^fe(x)) can be used to write the 
relations between metric distances fll9l) and (1201) in the following form 

r = S{t)Ij;\l/S), l = S{t)h{r/S), (22) 

where I^^ is the inverse function to Ik- At = we have r = I. 

The general formula for metric distance in the Friedmann model has the following 
form: 

Hq Jo n[z') 

where h{z) can be derived from the Friedmann equation in the following form: 



h{z) = ^p{z)no + {i-no){i + zy , (24) 

where = Ptot/Pcrit density parameter at the present epoch and p{z) = p/po is 

the normalized total density of all components. 

In case of the two-liquid dust + vacuum model (without interaction) the proper metric 
distance is given by the following formula: 



-no 



dz' 



and the external metric distance is equal to (for k = —1): 



(25) 



In the k = case 



In the k = +1 case: 



-h I r , 1 . (26) 



= , (27) 



K \ - ^ ^ 
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B. The "Angular Size Redshift" Relation 



The relation between metric distance l{z) and angular size has the following form: 

where d is the fixed size of the galaxy; x{z) =1{z)/Rhq and Rhq = cjH^^ Hq = 72 
(km/s)/Mpc. If is in arcsec and d in kpc then formula ( l29l) acquires the following form: 

e{z) = 0.0481 xdx (^j^) ■ (30) 




FIG. 4: Comparison of theoretical models with observational data for galaxies in the interval 
of absolute magnitudes M from -18 to -22 of the CiJ)viz_el-M18-22 sample. The dashed line 
corresponds to the theoretical model flm = 1-0 and = 0.0; the solid line — to the model 
= O.OandJ^y = 0.0; the dotted line — to Qm = 0.3 and = 0.7, and the dashed-and-dotted 
line — to = 0.0 and Oy = 1.0. 

For different combinations of Q^n and Qy nietric distance l{z) can be computed by 
formulae (I2S]-[2HD. 



C. Absolute Magnitudes 



To consider fixed intervals of galaxy luminosities, one must compute the absolute 
magnitudes of galaxies. In the Friedmann models the absolute magnitudes are computed 
using bolometric distance determined as Ihoi = l{z) x (1 + z). The absolute magnitude of 
the galaxy, Mj, in filter j can be computed from the following general formula: 

Mj = mj - 5 \og{l{z){l + z)} - 25 - AMj , (31) 

where distance l{z) is in Mpc and AMj = Aj + Kj + Ej is the correction to the absolute 
magnitude for extinction, redshift, and luminosity evolution. 



D. The "Surface Brightness-Redshift" Test 

The "surface brightness-redshift" test is a critical one, because, as Tolman f]\ was 
the first to point out that the J{z) relation is universal and the same for all Friedmann 
models. 

It follows from the definition of the surface brightness of an object that: 

Q2 - (1 + ^)4' ^'^^> 

where Jboi is the bolometric surface brightness; Fboi, the bolometric flux, and Jq, the sur- 
face brightness of the galaxy at z=0. In case of i-band we are dealing with the brightness 
observed in the given wavelength interval rather than with the bolometric brightness. 
Therefore below by Jboi we mean Jobs (observed in the the i filter). We then have: 

/i = /xo + 2.5 log(l + z)^ + Ki{z) + Ei{z) = fio + 2.5 log(l + z)", (33) 

where fx is the surface brightness measured in units mag"; Ki(z), the K-correction to the 
i-band surface brightness; Ei{z), the evolutionary correction to the i-band surface bright- 
ness; n = 4 + p = 4 + efc + ee, and p, the combined parameter of the surface-brightness 
evolution. 



V. EVOLUTION OF SIZES AND SURFACE BRIGHTNESS 
A. Parameter of the Evolution of Linear Sizes of Galaxies 

Figure H] compares the observational data for galaxies of the CiJjvizjel sample (here 
"el" indicates that the sample consists of elliptical galaxies) with theoretical models in 
the interval of absolute magnitudes M from -22 to -18. 

It is evident from Fig. H] that neither model passes through the median values for the 
given sample. This fact is usually explained by the evolution of linear sizes of galaxies 
and function /(z) is introduced: 



Q{z)obs = f{z) X Q{z)theor, 



(34) 



e, 



asec 
10 




FIG. 5: Finding the parameters of evolution for different models of galaxies with M from -20 to 
-18 for the CiJbviz-mod3-M 18 J2Q (a) and CiJ)viz-mod4:_M 18.20 (b) samples. The solid line in 
plots (a) and (b) corresponds to the models with = 1.0, Oy = 0.0 and ^Irn = 0.3, Oy = 0.7, 
respectively; the dashed line shows the fit for the entire sample, and the dotted line, the fit for 
median points. 



where Q{z)obs are the observed angular sizes and Q{z)theor are the theoretical angular sizes 
computed in terms of the model studied. The function usually has the form f{z) = (1+^)'^, 
where k is the parameter of evolution. 

We compute the parameter of the evolution of galaxy sizes in accordance with formula 
( EUj) . In Fig. 5 we present several plots for models with ^2^=1-0, 1^^=0.0 and ^2^=0. 3, 
f2y=0.7. We determine parameter k by applying the method of least squares both to the 
values for the entire sample and to median points. We perform our computations using 
mathematical package Microcal Origin 7.0. Note that the median is a statistically more 
stable parameter than the mean. As a result, we find the parameter of the evolution of 
galaxy sizes for four models both for the sample with M from -20 to -18 and for objects 
with M from -22 to -20. Table [T] lists the results obtained for angular sizes. 




FIG. 6: Evolution of surface brightness as a function of z for the CiJ)vizjmod3-M20J22 sam- 
ple. The dashed line corresponds to averaging over all points of the sample; the dotted line — to 
averaging over the median points; the solid line shows the theoretical evolution of surface bright- 
ness, and the vertical bars show the errors of the median values. 

Note that we set parameter d equal to 12 kpc for the sample of objects with M from 
-22 to -20 (compared to d=8 kpc adopted for galaxies with M from -20 to -18). This is 
due to the fact that bright galaxies (compared to galaxies with M from -20 to -18) have 
slightly larger sizes. That is why we slightly increased the "fixed galaxy size" . 



B. Surface Brightness Evolution 



The distribution of surface brightness values of the galaxies can be characterized by 
the J(0, z) profile for each z. In this paper we analyze the evolution of the surface 
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FIG. 7: Surface brightness evolution as a function of z for the CiJ)viz-mod3JM20J22 sample. 
The dashed line corresponds to averaging over all points of the sample, p = —0.8251; the dotted 
line corresponds to averaging over median points, p = —0.1144; the solid line shows the evolution 
of surface brightness, p = 0, and the vertical bars show the errors of the median values. 

brightness of galaxies, i.e., the central value of the J[z)@=o profile found using sExtractor 
and GALFIT codes. 

We determine the normalizing constant Jq by averaging the surface-brightness values 
for galaxies with z from to 0.5. The surface-brightness parameter n = 4 + p = 4+ k 
+ej is given by the following formula: 



where Jnorm is the normalized surface brightness: Jnorm = J/ Jo {Jo is the central surface 
brightness of the given galaxy) and p is the parameter of evolution. We determine param- 
eter n by applying the least-squares method both to all points and to the median values 
of normalized surface brightness. In Fig. [6] we compare the curves of the theoretical and 
observed evolution of surface brightness. 

For better visualization, these curves can be drawn in other axes. We use formula (l35i) 
to change to surface brightness : 



1 



(35) 



J, 



norm 





(36) 



whence it follows that 




(37) 



where Ax is the scale factor for converting pixels to acrseconds. In Fig. [7] we show the 
plots with different evolution parameters p in the (/x vs. lg(l + z)) axes, and /xq = 17.89 
mag" . 
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FIG. 8: Observed (the solid line) and modeled (the dashed line) binned (Az = 0.2) redshift 
distributions. 

The results of our analysis of the surface-brightness evolution of the galaxies 
of the HUDF field considered imply a surface-brightness evolution parameter of 
l? = — 0.83 ±0.1 (n = 4 +p = 3.17 ± 0.1) if computed using all points of the sample and 
p = —0.11 ± 0.1 (n = 3.89 ± 0.1) if computed by using only median points for galaxies in 
the interval of absolute magnitudes M from -18 to -20. 

VI. IDENTIFICATION OF POSSIBLE SUPERLARGE STRUCTURES IN 
THE RADIAL DISTRIBUTION OF HUDF GALAXIES 

The distribution of photometric redshifts of HUDF galaxies in the z- interval from 0.1 
to 6.5 can be searched for superlarge structures, which show up in the form of fluctuations 
of the number of galaxies within the corresponding bins of the distribution considered. 

A. Distribution of Photometric Galaxy Redshifts 

As a model distribution of galaxies to compare with the observed distribution, we use 
a uniform distribution of points inside a unit-radius sphere. We randomly assigned to 
each point an absolute magnitude in accordance with the Schechter luminosity function: 

= 0*(^)"exp d . (38) 

Our next step was to find the redshifts for modeled data points. To do this, we had to 
change from the radial distance unit to metric distance one. At /c = (a zero-curvature 



space) the outer and inner metric distances coincide, l{z) = r{z). We thus use the inverse 
relation z = l~^{z) and formula (!27|) to compute the redshift. To allow for selection due 
to the limit of telescope sensitivity, we limit the sample by apparent magnitude — it must 
not exceed 29"^. 

We fit the resulting model distribution by the following formula: 

dN = Ax''exp(-—Y dz, (39) 
V XqJ 

where the free parameters a, (3, and xq are inferred using the least squares method and 
A is the normalizing constant. 

TABLE I: Results of computation of the parameter k of the evolution of angular sizes of the galaxies 
for different cosmological models 



N 


Model 


Evolution parameter k 


M(mag.) 






based on all points 


based on median points 




I.l 


Qm = 0.0, ftv = 1-0 


-0.40 ±0.03 


-0.58 ±0.09 


-20-18 


1.2 


= i.o,nv = 0.0 


-0.99 ±0.04 


-1.09 ±0.05 


-20-18 


1.3 


Qn, = 0.3, J^v- = 0.7 


-0.79 ±0.03 


-0.91 ±0.05 


-20-18 


1.4 


Qni = 0.0, Qv = 0.0 


-0.67 ±0.03 


-0.89 ±0.18 


-20-18 


II.l 


Qm = 0.0, Qv = 1-0 


-0.49 ±0.06 


-0.60 ±0.08 


-22-20 


II.2 


Qm = 1.0, = 0.0 


-1.14 ±0.10 


-1.32 ±0.12 


-22-20 


II.3 


Qm = 0.3, Qv = 0.7 


-0.90 ±0.08 


-1.07 ±0.06 


-22-20 


II.4 


Qm = 0.0, Qv = 0.0 


-0.75 ±0.08 


-0.85 ±0.12 


-22-20 



B. Comparison of the Expected and Observed Distributions 

In Fig. 8 we present the observed and modeled (according to ( l39|) ) redshift distributions 
of HUDF galaxies. The parameters of the modeled distribution are a = 2.84, /3 = 0.48, 
xo = 0.015. 

We use the following quantity to measure the deviation of the observed number of 
galaxies from the theoretically expected numbers: 

^obs ^theor / An\ 

(Tabs = 77 , (40) 

' t/ieor 

where Ntheor is the expected (according to flHI?]) ) number of galaxies in the interval of 
redshift s from z to 2; + Az and Nobs is the number of galaxies observed in the same 
interval. 

The theoretically expected amplitude of fluctuations is characterized by the dispersion 
of Poisson noise, tXp = 1/ Ntheor, and the dispersion associated with correlated structures, 
which we computed by the following formula 1X3]: 
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z 

FIG. 9: Observed (the solid line), theoretical (the dashed line), and Poisson (the dotted line) 
deviations of photometric redshifts in Az = 0.2 bins. 

where 7 = 1.8; J2 = 72.0/[2^(3 - 7)(4 - 7)(6 - 7)]=!. 865; tq = 5 Mpc; r^ff = 
(3/47rr^Ar5') 3 is the effective radius corresponding to the volume of the interval; Ar 
corresponds to the dz layer, and S is the solid angle of the HUDF field. 

In Fig. 9 we present the plot of observed deviations aobs and theoretically expected 
deviations dp, dtheor from uniform distribution of galaxies. 

VII. DISCUSSION OF THE RESULTS AND THE MAIN CONCLUSIONS 

Our analysis of the classical cosmological test Q{z) shows that the choice of cosmo- 
logical models has a strong effect on the parameter of the evolution of linear sizes of 
galaxies. 

It is evident from Table [T] that the parameter of the evolution of linear sizes for galaxies 
with absolute magnitudes in interval from -20 to -18 varies from k= —0.40 ± 0.03 for the 
model with Q^n = 0.0 and Qy = 1-0 to k= —1.09 ± 0.06 for the model with = 1-0 
and Qy = 0.0. This parameter varies from —0.49 ± 0.06 to —1.32 ± 0.12 for galaxies 
with luminosities in the interval from -22 to -20. The inferred values of parameter k 
agree with the results of Bowens et al. [l^ for other samples of galaxies from HDF- 
S, HDF-N, GOODS, and HUDF. The above authors used a cosmological model with 
parameters Qm = 0.3 , fly = 0.7. For galaxies with luminosities in the interval of absolute 
magnitudes M from -22.38 to -21.07 the parameter of galaxy size evolution was found to be 
k= — 1.05±0.21. For our sample parameter k= — 1.07±0.06 in case of Qm = 0.3, Qy = 0.7. 
Note that the Q{z) diagram may become an efficient cosmological test when a reliable 
model of the evolution of galaxy sizes is developed. 

The surface brightness evolution parameter 

p = —0.11 ± 0.1 (n = 3.89 ± 0.1) inferred from the median points for galaxies in 
the absolute magnitude (M) interval from -18 to -20 requires further analysis. This is 
due to the fact that the K-correction to the surface brightness includes a combination of 




the K-correction to the flux and the K-correction to the angular size. 

An analysis of the distribution of HUDF galaxies reveals strong deviations of the ob- 
served number of galaxies from the number of galaxies expected for a uniform distribution. 
The observed irregularities correspond to a scale length of about 2000 Mpc. This may 
be due both to real superlarge structures and to hidden selection effects that show up in 
finding the photometric redshifts. This problem requires further analysis. 
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